The pseudogap phase in (TaSe4)2l 
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We have developed the mean-field theory of coexisting charge-density waves (CDW) and un- 
conventional charge-density waves (UCDW). The double phase transition manifests itself in the 
thermodynamic quantities and in the magnetic response, such as spin susceptibility and spin-lattice 
relaxation rate. Our theory applies to quasi-one dimensional (TaSe4)2l, where above the CDW tran- 
sition, thermal fluctuations die out rapidly, but robust pseudogap behaviour is still detected. We 
argue, that the fluctuations are suppressed due to UCDW, which partially gaps the Fermi surface, 
and causes non-Fermi-liquid (pseudogap) behaviour. 

PACS numbers; 71.45.Lr, 75.30.Fv, 75.20.En 

I. INTRODUCTION 

One dimensional interacting metals are unstable with respect to charge-density wave (CDW) formation. A strictly 
one-dimensional system cannot exhibit long-range CDW ordering above T = 0, although strong CDW fluctuations 
occur. Mean-field theories provide us with a finite Tc by neglecting the fiuctuations. However, increasing the dimen- 
sionality of the system by weak interchain coupling can lead to long-range ordered CDW at transition temperature 
Tc well below the mean-field transition temperature. Consequently, quasi-one dimensional systems were expected to 
possess CDW ground state, which was confirmed by experimentsi. 

Among the possible candidates, inorganic conductors as NbSea, TaSea, K0.3M0O3 and (TaSe4)2l have attracted 
much attention. Due to their quasi-one dimensional nature. CDW ordering and fiuctuation should be considered on 
equal footings, as was done in the pioneering work of Lee, Rice and Anderson^. Their results fitted a variety of results 
on the thermodynamic and magnetic properties of inorganic conductors in the fiuctuation regime slightly below and 
above T(~i4. However, the last member of the materials cited above, namely (TaSe4)2l, exhibits peculiar properties: 

1. Sharp increase of the static spin susceptibility above Tc in a wide temperature ranged. 

2. Smooth increase of the spin-lattice relaxation rate above T^. 

3. As opposed to conventional CDW, no sharp features in the spin lattice relaxation rate below T^. 

4. Thermal fluctuations die out rapidly with increasing temperature above T^. 
The first has been interpreted within the fiuctuating gap model^, whose application is questionable in light of the 
forth. The other two has not been addressed theoretically so far. 

Instead of using the fluctuating gap model, we propose an alternative model: the low temperature phase of (TaSe4)2l 
is dominated by CDW, but the response of the high temperature or pseudogap phase is attributed to unconventional 
charge-density waveai (UCDW). Unconventional charge density waves are density waves with a wavevector dependent 
gap, whose average vanishes on the Fermi surface ((A(k)) = 0). As a result, no modulation of the charge or spin 
density is expected in these systems, and the notion "hidden-order" applies naturally. The gap function has zeros 
at the Fermi surface, where low energy excitations, characterized by Dirac fermions, occur. The thermodynamic and 
transport properties of unconventional density waves have been worked out, for a review see Ref Q. Such phases 
were proposed for various materials, e.g. for the pseudogap phase in high superconductors^iS. 
^ , In the foUowings, we are going to discuss the coexistence of CDW and UCDW. Within the mean-field approximation, 
' we determine the phase diagram, and find the region of coexistence. Interestingly, a CDW gap can open in an existing 
UCDW phase, but the opposite cannot occur. The spin susceptibility and spin lattice relaxation rate keep on increasing 
when raising the temperature above the CDW Tc, in accordance with experimental data^. In addition, no critical 
divergence of the relaxation rate below Tc is found, as opposed to the usual CDW responsei^. Whence, the four basic 
properties of the pseudogap phase are satisfied by our CDW-I-UCDW model. For the first three, direct calculations 
are shown to evidence the agreement, while the last is argued to be obeyed since the number of normal particles is 
heavily suppressed above CDW due to the presence of UCDW condensate. To explore the features of our model in 
full glory, further studies are needed. 
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II. THERMODYNAMICS, PHASE DIAGRAM 

We consider the simple model Hamiltonian describing density waves given bjiii: 

H = Y^ [^(k)(aJ,^ak,<T - a+_(^^ak-Q.„) + A(k)aJ,,ak-Q,<T + A(k)a+_Q^^ak,<x (1) 

k,(7 

where a^^ and ak,(T are, respectively, the creation and annihilation operators of an electron of momentum k and 
spin a. In a sum with prime A:^. runs from to 2kp (kp is the Fermi wavenumber), Q = (2fc^, 7r/6, tt/c) is the best 
nesting vector. Our system is based on an orthogonal lattice, with lattice constants a, 6, c toward direction x,y,z. 
The system is anisotropic, the quasi-onc-dimcnsional direction is the x axis. The linearized kinctic-cncrgy spectrum 
of the Hamiltonian is: 

^(k) = VF{kx — kp) — 2tbCos{kyb) — 2tcC0s{kzc) — fi. (2) 
A(k) is the density wave order parameter, and is determined from the self-consistency condition: 

A(l)-^P(k,l)(a+,ak-Q,.>, (3) 

k,(T 

where P(k, 1) is a linear combination of the interaction matrix elements, in which we considered at most nearest 
neighbour interactionsii. Since we are focusing on the coexistence of CDW+UCDW, we take the kernel of the form 
P(k, 1) = —Pq — Pis'm{bky)sm{bly), with both Pq, Pi > 0, using the analysis of the electron-electron interaction in 
quasi-one dimensional systemsii. The sign of the matrix elements assures charge-density wave formation. Had we 
chosen a gap with cosine function or with the z component of the wavevector would not alter our results. Based on 
this kernel, the gap is given by 

A(k) = Ao + Aie"^ sm{bky). (4) 

Here we have skipped the overall phase of the condensate since it turns out to be irrelevant for most of our discussion. 
The value of the relative phase of the gap coefficients, should be determined by minimizing the free energy of the 
system with respect to it. It turns out, that coexisting solution can occur only for = or tt/2 (by adding integerXTr 
describes the same physics). Among these, the former solution has always higher energy than the pure phases, namely 
CDW with Ao 7^ and Ai = or UCDW with Aq = and Ai ^ 0. The latter provides us with energetically stable 
coexistence, and the gap now is of the form: A(k) = Aq + lAi sm(bky). The coupled gap-equations are obtained as 

7r/2 W 


tt/2 W 

„ ^ f 2A{bky) f ,^ . 2„, , tanh (/3£;/2) 
Ai^po^^iAi J -^-Jil J d^sm\bky) 



where 

^^(k) = ^e{^) + Al + Ajsm^{bky) (7) 

is the energy spectrum, W = vpkp is half of the bandwidth, po is the density of states in the normal state per spin 
at the Fermi energy. This spectrum implies a finite energy gap as long as Aq remains finite. The phase diagram 
is depicted in Fig. ^ Coexistence region is found only around 2Pq < Pi, because this implies that the transition 
temperature to CDW or UCDW is close to each other, hence the two phases can compete. The gap-equations have 
been solved numerically and the result is shown in Fig. [3 The system at low temperatures can enter into a CDW 
or UCDW or CDW+UCDW phase. In the following we are focusing on the latter region of coexistence. By lowering 
the temperature, first the UCDW gap opens at the transition temperature Td = 27T4^exp(— 4/Pipo)/7r given by the 
weak-coupling solutionii of Eqs. I5I6I 7 = 1.781 is the Euler constant. It follows the T dependence of the pure solution 
until the CDW gap opens at Tco- Then Ai slightly decreases and saturates to a finite value while Aq grows in a 
mean-field manner with decreasing T. The opposite scenario, the opening of UCDW gap upon an existing CDW one 
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FIG. 1: (Color online) The mean-field phase diagram is drawn for UCDW coupling Pipo ~ 1 with transition temperature Td, 
Tci/W ~ 0.0202. The dashed line stands for the pure CDW transition in the absence of UCDW. The transition along solid 
lines is of second order. Similar phase diagram is obtained for other values of Pi. 
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FIG. 2: (Color online) The temperature dependence of the gaps are plotted for various interaction strength: Pipo = 1 and 
Popo = 0.49 (blue solid line), 0.48 (red dashed line), 0.47 (black dashed-dotted line) and 0.46 (magenta dotted line). The gaps 
are normalized by the gap maximum of a pure UCDW at T = 0. The curves starting at the highest Tc belong to Ai. 



is excluded within the present model, because no such solution of Eqs. Q-® was found. At T = 0, the coexistence 
region is located to 4/(8 + i'LPo) < Po/Pi < 1/2- The larger Pi, the wider this region is. In the coexistence region, the 
amplitude of the charge-density oscillations is determined solely by Aq, while in the pure UCDW phase, no modulation 
is found. 

Since within the mean-field theory, we are dealing with non-interacting quasiparticles, their entropy can be calculated 
from the knowledge of the energy spectrum and yields to 

SiT) = -2^ (/ln(/) + (1 - /)ln(l - /)) , (8) 

k,(T 

where / ~ l/(exp(/3i?(k)) + l) is the Fermi function. By taking its derivative with respect to temperature, the specific 
heat can be evaluated. This double phase transition is manifested in the specific heat as well, which exhibits two 
peaks at the transitions to UCDW and CDW, as can be seen in Fig. O 



FIG. 3: (Color online) The temperature dependence of the specific heat is shown for Pipo — 1 and Popo ~ 0.49 (blue solid 
line), 0.48 (red dashed line), 0.47 (black dashed-dotted line) and 0.46 (magenta dotted line). The curves are normalized by the 
the normal state T linear specific heat, 7s is the Sommerfeld coefficient. 



III. MAGNETIC RESPONSE 

The magnetic response of the model studied in the foUowings includes the long wavelength static spin susceptibility 
and nuclear spin-lattice relaxation rate. These are related to the single particle density of states, hence we begin with 
the exploration of its properties. We find 



tt/2 



per spin direction. Here Ncdw{x) = \x\/Vx'^ ~ lQ{\x\ — 1) is the density of states in a conventional DW, Nucdw = 
8(1 — \x\)2xK{x)/'K + Q{\x\ — 1)2K{1/x)/tt is the density of states in UDW. Whence, the total density of states 
in the coexistence region is characterized by a finite energy gap Ag stemming from the appearance of CDW, and a 
logarithmic singularity at ±-\/A| + Af reminiscent to that in UCDW. As soon as Ai disappears, this turns into a 
square-root divergence. These features can be observed in Fig. ^ 

From this, the spin susceptibility and the spin-lattice relaxation rate are calculated as 

X(r)_ 7mE)df^^^ ^^^^ 



dE, (11) 



Xo J Po dE 

— OO 
OO 

R{T) f N^{E) df 



Ro J pI dE 



where xo and Rq are the respective metallic state values. The first expression applies to various kinds of condensates. 
The latter one is determined from the short wavelength spin response (q ~ Q), because dominant contribution comes 
from wavevectors close to the ordering (nesting) vector. This explainsifiji^ the absence of coherence factors even in 
the presence of a constant gap as opposed to the type II coherence factors in superconductors^^. By inserting the 
solution of the gap equation into Eqs. (|10|I - H11() . we obtain the temperature dependence of these quantities. At low 
temperatures, activated behaviour is seen because of the finite CDW gap. Upon increasing T, a change of slope is 
exhibited at Tco, and then the curves follow that in a pure UCDW. In the spin-lattice relaxation rate, a small peak 
is observable right below Td- In Maniv's original expressioni^ for i?(T) in a conventional DW, the magnitude of the 
peak was controlled by Aln(8A/wL), ujl being the Larmor frequency. In our case, the divergence in the density of 
states is logarithmic, hence the Larmor frequency can be taken to zero. Consequently, no sharp features are observed 
close to Tco as seen in Fig. |S1 Nevertheless, a small peak remains right below Td as a reminder to the density wave 
nature of the condensate. 
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FIG. 4: (Color online) The temperature and frequency dependence of the quasiparticle density of states is shown for Pipo = 1 
and Popo ~ 0.48. The temperature dependence enters through Ai{T). 
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FIG. 5: (Color online) The temperature dependence of the spin susceptibility is shown in the left panel for Pipo = 1 and 
Popo = 0.49 (blue solid line), 0.48 (red dashed line), 0.47 (black dashed-dotted line) and 0.46 (magenta dotted line). The right 
panel shows the spin- lattice relaxation rate. Here the deviation for various Po's is smaller than for the spin susceptibility. 

IV. APPLICATION TO (TaSe4)2l 

In mean field theories, attractive interaction on lattices leads to CDW formation in low dimensional systems. 
In a strictly onc-dimensional chain, however, CDW fluctuations eliminate any long range order above T = 0. In 
reality, finite interchain coupling increases dimensionality and restores long range order at a temperature, which is an 
appreciable fraction of the mean-field value. 

This is why quantum fluctuations are very important in studying low dimensional systems such as K0.3M0O3, TaSa 
and other CDW materials. Lee, Rice and Anderson^ investigated the effects of fluctuations on the Peierls transition 
in one dimension by taking a functional integral over variations in the order parameter. The success of their theory 
became apparent when applied to the aforementioned CDW materials. They found the suppression of the mean- 
field Tc by the fluctuation. Moreover, the fluctuating region extended to high above the transition temperatureii. 
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and generated pseudogap behaviour. Its manifestation includes the smooth increase of the spin susccptibihty before 
reaching the PauU hmit^. 

As opposed to this, (TaSe4)2l exhibits special phenomena. When raising the temperature above the charge-density 
wave transition temperature at T = 265 K, thermal fluctuations disappear in X-ray measurements^, but the character- 
istic pseudogap behaviour of increasing spin susceptibility^ and spin-lattice relaxation rate^ remains intact. Whence 
these cannot be related to fluctuations. In an effort to resolve this contradiction, we propose an alternative sce- 
nario: the material does not enter into its normal metallic state above 265 K, but remains in UCDW phase. This 
is characterized by a gap, which vanishes on certain subsets of the Fermi surface. At other points, quasiparticles 
experience the full energy gap. Fluctuations are still present, but their effect is reduced by the available small portion 
of non-condensed particles ( 50 — 95% of the original electrons may belong to the UCDW condensate, as seen in Fig. 
El through the relation: condensate density^ 1 — x(^)/Xo)- In this respect, our model also accounts for the reduced 
electronic density of states at the Fermi energjii^. Pseudogap behaviour is produced by excitations living around 
the gap nodes. In the preceding sections, we have successfully demonstrated that CDW-I-UCDW theory reproduces 
the basic features seen in magnetic probes. The slope of both the spin susceptibility and the spin lattice relaxation 
rate changes when passing through the CDW transition, and both keep on increasing with increasing temperature. 
Moreover, a critical divergence of the spin lattice relaxation rataiS at the CDW transition temperature occurs usually. 
As opposed to this, no sharp feature in R{T) has been detected experimentally^ close to Tco, in accordance with our 
findings, as seen in Fig. [S] In Fig. El we plot the spin-lattice relaxation rate as a function of x^iT). The deviation 
from the Korringa relation, R{T)/x^{T) ~ const is similar to that found in (TaSe4)2l^- 




FIG. 6: (Color online) The Korringa plot of the relaxation rate for Pipo — 1 and Popo = 0.49 (blue solid line), 0.48 (red dashed 
line), 0.47 (black dashed-dotted line) and 0.46 (magenta dotted line). The straight green dashed line is R{T)/Ro = (x(2^)/xo)^- 



We speculate, that the transition from UCDW to normal metallic state takes place above 430 K, when (TaSe4)2l 
starts to decompose. As to the conventional CDW transition, it is surely overwhelmed by the large phononic contri- 
bution in specific heat, although some anomalies are reportedi^. 

As seen in Fig. Q] the coexistence is localized into a narrow region between pure CDW and UCDW. When the 
system is forced to leave this regime, it should enter into either a pure CDW or UCDW. By the application of 
chemical or mechanical pressure (which, by increasing the interchain hopping, destabilizes DW phases) or by the 
onset of external magnetic field, such a transition can be induced. The UCDW transition can in this way be brought 
to the temperature range of feasible experiments. 

Within our weak-coupling mean-field analysis, agreement between theory and experiment is achieved only on 
qualitative levels. The experimental value of the single particle CDW gap is^ around 1200 — 1500 K, which implies 
the gap over transition temperature ratio to be 5 — 6, much larger than the weak-coupling BCS value 1.76. We note, 
that taking imperfect nesting into accountiii^ would significantly enhance this ratioiS, and quantitative agreement 
could be reached. This is beyond the scope of the present paper, but further research in this direction is desirable. 
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V. CONCLUDING REMARKS 



We have studied the coexistence of CDW and UCDW within the mean-field approximation. Stable coexistence is 
found only for complex order parameter, when the relative phase of the CDW and UCDW gaps is ±7r/2. By lowering 
the temperature from the normal state, first UCDW opens at Td. Than CDW shows up at T^o upon the existing 
UCDW phase. The CDW (UCDW) gap increases (decreases) monotonically with decreasing temperature until T ~ 0. 
Periodic oscillation of the charge density is expected only in the CDW+UCDW phase. The opposite scenario, the 
opening of UCDW gap in a CDW phase has not been found in our model. The density of states possesses a clean 
gap as long as the CDW phase is present, and the square root singularity characteristic to CDW is replaced by the 
weaker logarithmic one of UCDW. The specific heat exhibits two jumps at the corresponding transition points. The 
spin susceptibility and spin-lattice relaxation rate increases steadily without any signs of saturation or coherence peak 
through Tco because of the logarithmic nature of the quasiparticle peak in the density of states. 

This study is, however, not just of academic interest. Salient features measured in the pscudogap phase of (TaSc4)2l 
can be explained with our CDW-I-UCDW theory. These are: the increase of both the spin susceptibility and spin- 
lattice relaxation rate with temperature above T^o, the featureless behaviour of the latter right below Tco, and the 
rapidly decaying thermal fluctuations around the CDW transition. In this manuscript, we have demonstrated that 
these phenomena are readily described by our coexisting CDW and UCDW model. Further studied in this direction 
are highly requested in order to reach decisive conclusions. 
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